Empirical formulas for the fermion spectra and Yukawa matrices 
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We present empirical relations that connect the dimensionless ratios of fermion masses for the 
charged lepton, up-type quark and down-type quark sectors: \V US \ ~ [~1 5 ~ V^T\ * ~ 3 
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Explaining these relations from first prin- 



and || |^-| 

ciples imposes strong constraints on the search for the theory of flavor. We present a simple set 
of normalized Yukawa matrices, with only two real parameters and one complex phase, which ac- 
counts with precision for these mass relations and for the CKM matrix elements and also suggests 
a simpler parametrization of the CKM matrix. The proposed Yukawa matrices accommodate the 
measured CP-violation, giving a particular relation between standard model CP-violating phases, 
(3 = Arg [2 — e~ !7 l . According to this relation the measured value of /3 is close to the maximum 
value that can be reached, /3 max = 30° for 7 = 60°. Finally, the particular mass relations between 
the quark and charged lepton sectors find their simplest explanation in the context of grand unified 
models through the use of the Georgi-Jarlskog factor. 



I. INTRODUCTION 

Any theory of flavor must explain the fermion mass 
hierarchies as well as the quark mixing angles. Unfortu- 
nately, few patterns have been found in the measured val- 
ues of fermion masses and mixing angles that can guide 
us in the search for an underlying theory of flavor. One 
of these, which has been known since 1968 [lj, is the well 
known empirical relation between the down-quark mass, 
the strange-quark mass and the Cabibbo angle, 
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(1) 



This relation has driven the development of theories of 
flavor over more than three decades, starting in 1977 with 
the first attempt to explain it using family symmetries 
0, L| • Another quark mass relation that has been known 
for some time is, 



(2) 



Inspired by these two relations, many of the theories 
of flavor proposed to date have focused on generating 
Yukawa matrices that are polynomial in powers of A, 
A ~ \V US \, with coefficients of order 1 Q. There is a 
third famous relation. It was argued as early as in 1979 
that at momenta larger than 10 15 GeV quark and charged 
lepton masses are related by, 
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m d = 3m e 



to m = 3m s , 



(3) 



An ingenious method was proposed to account for this 
relation by the use of SU(5) Clebsch-Gordan coeffi- 
cients Q- Other than these relations, it is usually claimed 
that the fermion masses follow scaling laws of the form 
(md,m s ) ~ (A 4 ,A 2 )mf, in the down-type quark sector, 
(m u ,m c ) — (A 8 ,A 4 )m t in the up-type quark sector and 
(m e ,TO M ) — (A 4 ,A 2 )m T in the charged lepton sector. As 
can be easily checked, however, these scaling laws are 
qualitative and do not survive a precision analysis. 

The measurement of the top quark mass in 1995 and 
the continuous improvement in the extraction of other 
quark masses during the last decade motivate a more sys- 
tematic search for precise empirical relations between di- 
mensionless ratios of fermion masses in each fermion sec- 
tor. There are six independent fermion mass ratios of this 
kind, two for each fermion sector. It is possible for hidden 
regularities to manifest themselves more clearly through 
higher order dimensionless ratios of fermion masses, i.e. 
ratios of the form ra 2 ,/ '(mf,m c ) or m a /(mlm c ). Indeed, as 
we show in this paper, there are some interesting patterns 
underneath the measured values of the fermion masses. 
These new relations, which are not merely qualitative, 
are the following, 
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(5) 



We expect these two basic parameters, which we denote 
hereafter by A ~ \V US \ and 9 ss | \V c b\ / \V US \, to be con- 
nected with the fundamental parameters of the underly- 
ing theory of flavor. 

This paper is organized as follows. We begin in Sec.lTTl 
by systematically searching for correlations between di- 
mensionless mass ratios in different fermion sectors up to 



order 3, i.e. up to ratios of the form m? a / '(m%m c ). We re- 
view in the appendix the calculation of lepton and quark 
running masses which are used in Sec. [n] In Scc lIIII we 
analyze Yukawa renormalization corrections that affect 
the studied mass relations when evolved with the renor- 
malization scale, especially to the ratios including third 
generation fermion masses. In Sec. II VI we show that, as a 
consequence of these new empirical formulas the fermion 
mass hierarchies can be expressed as a function of two 
basic parameters, A and 9. In Sec. |V] we show, neglect- 
ing CP-violation, that the absolute values of the CKM 
mixing matrix elements can also be expressed as simple 
functions of the basic parameters A and 9. In Sec. I VII we 
propose, neglecting CP-violation, a simple reconstruction 
of the quark Yukawa matrices that accounts for the corre- 
lations found in the previous sections. In Sec. lVIII we in- 
troduce CP-violation in the textures proposed in Sec. I VII 
and study its predictions for the CP-violating parame- 
ters. In Sec. IVIlil we study the precision predictions for 
the lighter quark masses, CKM elements and charged lep- 
ton masses arising from the texture proposed in Sec. lVIl 
In Sec. lIXI we point out that the simplest solution to ac- 
count for the relations between the charged lepton sector 
and the quark sector can be found in the the extension of 
the Standard Model SU(3) C X SU(2) L x U(1) Y symme- 
try to the SU(5) symmetry of Georgi and Glashow. In 
Sec. |X] we speculate about the characteristics of under- 
lying flavor models that can reproduce these empirical 
mass relations. 
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FIG. 1: The fermion mass spectra. 



of the three fermion sectors: up-type quark, down-type 
quark and charged lepton to be expressed as a simple 
polynomial functions of the flavor parameters: A, 9, ■ ■ •, 



II. CORRELATIONS BETWEEN 
DIMENSIONLESS FERMION MASS RATIOS 

In this section we will look for patterns in the dimen- 
sionless mass ratios of running fermion masses. Other 
than the fact that the first fermion generation is lighter 
than the second and this is lighter than the third genera- 
tion, there are no other evident regularities in the fermion 
mass spectra, as can be observed in Fig. ^ Based on 
the experimental fact that the third generation is much 
heavier than the first and second generations and that 
the quark mixing angles are small we hope that there is 
a simple mechanism of flavor breaking which generates at 
some higher energy scale a simple structure in the nor- 
malized Yukawa matrices. If this is the case it is plausible 
that at such scale the normalized Yukawa matrices have 
the form, 
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(6) 



where A, 9, ■ ■ ■ represent generically some perturbative 
flavor breaking parameters, i.e. A, 9, ■ ■ ■ <C 1, directly 
related to the underlying theory of flavor. We note that 
in many flavor models proposed in the literature the fla- 
vor breaking is parametrized by a unique parameter A. 
Therefore we expect the fermion mass ratios in each one 
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Let us assume to simplify the discussion that there are 
only two flavor breaking parameters: A and 9. In this 
case it would be possible to solve the previous system of 
equations and obtain expressions for A and 9 as a function 
of the fermion mass ratios, 
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This can be done for each fermion sector separately. This 
makes plausible that underlying patterns manifest more 
clearly in higher order mass ratios, even though these 
can be expressed as a function of the six basic fermion 
mass ratios. When searching for mass relations between 
different fermion sectors, it is convenient to calculate ra- 
tios of running fermion masses at a common renormal- 
ization scale. If there are regularities in the underlying 
Yukawa matrices, these will be manifested more clearly 
in the ratios of running fermion masses, not in the ratios 
of physical masses. Using the running masses that we 
have calculated in the appendix we obtain dimensionless 
mass ratios in the charged lepton, up-type quark and 
down-type quark sectors. We calculate ratios of order 
1, cW, order 2, c} 2 \ and order 3, c^. These ratios are 



generically of the form, 
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where a,b,c — 1, 2, 3 are generation indices. Our numer- 
ical results are shown in table [I] We have also included 
uncertainties for the mass ratios, Ac, calculated using 
Ac = \dc/dm a \ Am a , where Am a are the uncertainties 
in the determination of running fermion masses. The 
measured quark and charged lepton masses used as an 
input in our calculations are explained in detail in the 
appendix. We have compared analogous c coefficients in 
the three different fermion sectors, looking for simple cor- 
relations of the form, c p — rc q or c p = c 7 q where r is a low 
integer number and p,q — l,u,d denote similar ratios in 
the charged lepton, up-type or down-type quark sectors. 
We have first searched for correlations between the 
mass ratios in the up-type quark and down-type quark 
sectors. We have only found two clear correlations. The 
first correlation appears for the order one coefficient in 
entry I of table HJ The correlation appears between the 
ratios, 
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0.211 ±0.033, 



0.225 ±0.018, 



(12) 
(13) 



and the Cabibbo angle \V US \. These ratios have uncer- 
tainties respectively of the order ±16% and ±8% of the 
central values. It is convenient to show this correlation 
in an alternative form, which makes it more manifest, 
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1.06 ±0.25 



(14) 



This correlation has been known for some time. Curi- 
ously, we also find an interesting correlation with the 
analogous ratio in the charged lepton sector, 
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3.06 ±0.4 



(15) 



The ±15% uncertainity in the calculation of this ratio 
comes from the uncertainty in the determination of the 
lighter quark masses. This indicates that the following 
ratio in the charged lepton sector, 
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= 0.20648 ±0.000002 



(16) 



gives a numerical value very close to the Cabibbo angle 
and to the ratios in Eas.lrHandlT^l It was first pointed 
out in 1979 by Georgi and Jarlskog [5j that at momenta 
larger than 10 15 GeV quark and charged lepton masses 
seem to be related by, 



m d = 3m e 



to m = 3m s , 



(17) 



and that this relation could be explained by the use of 
SU(5) Clebsch-Gordan coefficients. We want to empha- 
size that the previous correlations indicates that indeed 
there is a very precise relation between the Cabibbo an- 
gle and the ratios of the fermion masses of the first and 
second generation, 
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We will see in Sec. lIIII that. as a very good approximation, 
this relation is renormalization scale independent. We 
have found only one more simple correlation amongst the 
dimensionless mass ratios shown in table [I] This appears 
for the order three coefficient shown in the entry XI in 
table |U It is convenient to take the square root of the 
numbers shown in the table. In the up-type and down- 
type quark sectors we obtain, 




= 0.073 ±0.023, 



= 0.114 ±0.039. 



(19) 
(20) 



Both ratios have an an important uncertainty, approxi- 
mately ±32% of the central value, coming from the un- 
certainties in the extractions of the lighter quark masses. 
It is convenient to quantify the correlation by taking the 
ratio, 
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1.5 ±1.0. 



(21) 



We note that there are two integer numbers, {1,2}, inside 
the error bars which could give us a simple correlation. 
Furthermore, using in Eg. 1201 a recent lattice extraction of 
the strange quark mass mentioned in the appendix J20j 
instead of the sum rules extraction the ratio in Eq. I2T1 
turns out to be closer to 1 while the uncertainity is re- 
duced, 



~(m\ at ) 3 ~ 
m\m d _ 


2 


r sn 
m c 


m^m u 



= 0.8 ±0.5. 



(22) 



This lattice extraction of the strange quark mass has 
not been used throughout the main text because it has 
not yet been confirmed by other lattice QCD collabora- 
tions. If we compare with the analogous ratio between 
the down-type sector and the charged lepton sector we 
obtain, 
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= 7.5 ±2.6 



(23) 



The uncertainty is again important, approx. ±34% of the 
central value, but it clearly points out that there are four 
integer numbers, {6,7,8,9}, inside the error bars which 





Charged leptons 


Down-type quarks 


Up-type quarks 


I 


rrie/m^ 


(4.73711 ± 0.00007) x 10" 3 


m d /m a 


(4.4 ±1.4) x 10" 2 


m u /m c 


(2.6 ±0.8) x 10 -3 


II 


rn^/rnr 


(5.882 ±0.001) x 10~ 2 


m a /m b 


(2.4 ±0.4) x 10~ 2 


m c /m t 


(3.7 ±0.6) x 10~ 3 


III 


m e m^/m 2 


(1.6390 ±0.0006) x 10~ 5 


m d m 3 /m 2 b 


(2.5 ±1.0) x 10" 5 


m u m c /rn 2 


(3.65 ±1.5) x 10 -8 


IV 


mi/m^mr 


(1.3199 ±0.0003) x 10~ 6 


m d /m a m b 


(4.7 ±2.5) x 10~ 5 


ml/mcmt 


(2.52 ±1.4) x 10" 8 


V 


ml/m e m T 


12.417 ±0.002 


m 2 Jm d m b 


0.53 ±0.26 


m 2 /m u rn t 


(1.45 ±0.71) 


VI 


m%,/m 2 ^m T 


(6.253 ±0.001) x 10~ 9 


m d /m 2 m b 


(2.1 ±1.8) x 10" 6 


m^/m^mt 


(6.5 ±5.8) x 10" 11 


VII 


3 / 2 

m e /m^m T 


(3.678 ±0.001) x lO" 10 


m d /m s ml 


(5.0 ±3.7) x 10" 8 


m^/m c m 2 


(2.44 ±2.0) x 10~ 13 


VIII 


ml/m 2 m T 


2621.2 ±0.5 


ml/m d m b 


12 ±10 


m^./m 2 i mt 


(560 ± 455) 


IX 


m e ml/ml 


(9.640 ± 0.005) x 10~ 7 


m d m 2 s /ml 


(6.0 ±3.5) x 10~ 7 


m u m 2 /rn'l 


(1.4 ±0.8) x 10~ 10 


X 


mim^/mr 


(4.567 ± 0.002) x 10~ 9 


m 2 d m s /ml 


(2.7 ±1.6) x 10~ 8 


m 2 l m c /m'l 


(3.5 ±2.4) x 10 -13 


XI 


ml/m l ,m T 


0.7304 ± 0.0002 


m z s /m d ml 


(1.3 ±0.9) x 10" 2 


m^,/m u m t 


(5.4 ±3.5) x 10" 3 



TABLE I: Dimensionless fermion mass ratios in the charged lepton, up and down-type quark sectors calculated from measured 
values as explained in the text 



could give us a simple correlation. It is interesting to 
check the values predicted by multiplying by the inverse 
of these integer factors the coefficient in the charged lep- 
ton sector. For instance, multiplying by 1/8, 1/9 and 
1/10 we obtain the following values, 



0.10681 ±0.00001, (24) 
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= 9 

1 

To 
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ml 


m e m 2 



0.09495 ±0.00001, (25) 



0.08545 ±0.00001. (26) 



If we compare these values with the values of the coef- 
ficients in Eqs. ^5] and |2U] we observe that the integer 
factors 9 and 10 give us a number which is compatible 
with the error bars of the coefficients in the up and down- 
type sector simultaneously. We find specially interesting 
the appearance of the factor 9 because as we will see in 
section ITXl there is already a simple solution to explain 
this factor, the Georgi-Jarlskog factor in GUT theories. 
These results indicate that there may be a second pre- 
cise relation between the low energy ratios of the fermion 
masses of the first, second and third generations, 
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(27) 



The present uncertainities in the extraction of the lighter 
quark masses do not allow us to determine if the relation 
works at the 1% level or just at the 40% level. Surpris- 
ingly we must emphasize that the exact empirical relation 
as given by Eq. |^| is inside the lcr experimental uncer- 
tainities for all of the three coefficients as calculated in 
Eas. lT^l[2l7l andl2*5l These can be observed more clearly in 



Fig. [3 We hope that near future improvements in the ex- 
traction of the ligther quark masses, by the use of lattice 
QCD methods, could test with precision this empirical 
formula. We note that there are only six independent 
mass ratios. We have already found two simple and pre- 
cise correlations linking the six of them. Therefore there 
cannot appear new correlations for other dimensionless 
mass ratios that cannot be expressed as a function of 
these two. 



III. FERMION MASS RATIOS AND THE 
YUKAWA SCALE 

In Sec. Hill we searched for correlations between di- 
mensionless ratios of fermion masses at low energies. 
Nonetheless, it is known that the top quark Yukawa cou- 
pling is of order 1 , since the top mass is of the same or- 
der than the electroweak scale. This implies that Yukawa 
coupling corrections cannot be ignored in the renormal- 
ization of the third generation fermion masses to very 
high energies. In models where the Higgs fields are flavor 
independent approximate solutions that relate mass ra- 
tios at different renormalization scales fi and (1q are given 
by 
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(28) 
(29) 
(30) 



In the case of the SM and the MSSM these effects can be 
calculated approximately from available one-loop renor- 
malization group equations 6] . For the Standard Model 




10 " 10 " 10 



Q (GeV) 

FIG. 2: Renormalization scale evolution of the coefficients 
8 t , 9b and 6^' 9 ' 10 couplings and their uncertainities according 
to the SM RGE equations. The three coefficients and their 
uncertainities shown in the plot were calculated using Eqs. \1!A 
\2(A and \2JJk£b\ The ligther quark masses used were extracted 
using sum rules. 



we obtain £j = £ t = £ and £ is denned by, 




£«Exp 



since mb/rrit ~ 2 x 10 2 this is approximately, 



t»[JL 



(31) 



(32) 



Moreover the tau lepton Yukawa renormalization factor 
is very small, 



£r«£ 



fe) 2 



(32^2 10 4 ) 



(33) 



We arc interested in evaluating how the correlations 
found in Sec. ^ evolve with the renormalization scale 
when including Yukawa corrections. To this end let us 
define the dimensionless ratios, 



n = 



m\m d 



(34) 
Their evolution, using Eas. 1281 1291 and 1301 is given by, 



e 



e (35) 



If we assume that fi/fio = Mg/Mz ~ 10 14 we obtain 
£ w 1.36. If we extrapolate the mass relations up to the 
GUT scale using SM RGEs we obtain, 



Mo 



and 



2.8 ±1.0, 



3.5 ±2.6, 



(36) 



(37) 



which must be compared with the low energy ratios cal- 
culated in Sec. Ull 



f =1.5 ±1.0, 

y * / M Z 



and 



7.5 ±2.6. 



(38) 



(39) 



M z 



Therefore the renormalization up to the GUT scale seems 
to spoil the mass correlation between the up and down- 
type quark sector. These results are summarized in Fig. [21 
and they may indicate that the Yukawa scale, the scale 
where the Yukawa couplings are generated, is an inter- 
mediate scale much lower than the GUT scale or alterna- 
tively that is not correct to use SM RGEs in the evolution 
of the fermion masses up to the GUT scale. If we assume 
that the couplings evolve according to the MSSM RGEs 
the results depend on tan/3, the ratio of Higgs expecta- 
tion values in the MSSM. We obtain, 



6 « e t /3 « Exp 



& ~ & ~ Exp 



1 Tnfii 
Mo 



16tt 2 



, tp * 1, (40) 
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4tt 2 Vmo 



, tfi > 1. (41) 



Therefore we obtain the following scaling factors for low 
tan/3, 






4/3 



-2/3 



(42) 
(43) 



Here £ was defined in Ea.l3*2land using /x/^o = Mg/Mz 
we obtain £ 4 / 3 rs 1.5 and £ -2 / 3 w 0.81. For large tp we 
obtain, 



^b_ \ ,_ 8 /3 
Mo 



(44) 
(45) 



Here £" 8 / 3 w 0.44 for /z//z = M G /M Z . Therefore in 
the MSSM, for both cases low and large tan/3, we can- 
not extrapolate the mass relations to scales as high as 



the GUT scale without spoiling the successful low energy 
mass relations. This again may be an indication that the 
Yukawa scale is not so far from the electroweak scale. We 
would like to point out that the ratios of first to second 
generation fermion masses; md/rUg, m u /m c and m e /m^, 
receive tiny Yukawa renormalization factors. Therefore 
the mass relation, 
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= 3 


m e 
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1/2 



(46) 



can be considered renormalization scale independent. To 
sum up, this analysis indicates that the second empirical 
relation, as given by Eg. 1271 may be optimal at some in- 
termediate or low energy scale. Nevertheless, the present 
uncertainities in the ligther quark masses are not small 
enough to allow us the determination of the scale at which 
this empirical formula is optimal. 



IV. THE FERMION MASS HIERARCHIES 

The empirical formulas found in the previous section 
can be simply understood if the fermion mass hierarchies 
are expressed as a function of two real parameters that we 
will denote hereafter by 9 and A. Let us assume that the 
ratios of running masses can be written in the following 
form, 



(m d ,m s ) = 
(m u ,m c ) = 

(m e ,m fi ) = 



(9\ 3 ,9\)m b , 

(ex 6 ,ex 2 )m t , 



-6»A d , 3 0A | ni., 



(47) 
(48) 

(49) 



We can easily prove that if this is the case we obtain 
immediately the correct empirical mass relations, 
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-,1/4 
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1/2 



in) 



= 3 

1/2 



1/2 



(50) 



mim v 



1/2 



(51) 



In other words we can say the the hierarchies in Eg s. 1471 - 
IH] solve the Eqs. I5UH5T1 We pointed out before that A 
corresponds approximately to the Cabibbo angle. The 
parameter 9 may be considered a new flavor parameter 
that seems to suppress, in all three fermion sectors, both 
first and second generation masses in the same amount 
with respect to the third generation. The fact that 9 and 
A connect different fermion sectors and that the correla- 
tions we have found work at a quantitative level suggest 
that 9 and A could be directly related to the underlying 
theory of flavor. We expect that any theory of flavor must 
be able to explain these correlations from first principles, 
and perhaps provide a prediction for 9 and A. 

Although the important uncertainties in the current 
extraction of lighter quark masses make impossible to 



know to what extent the relations in Eos.l50landl51lhold. 
it is plausible to study the implications derived of assum- 
ing that these relations are exact or almost exact. If this 
were the case we are lead to assume that the charged 
lepton sector is providing us the most precise determina- 
tion of the basic flavor parameters A and 9, A « 0.21 and 
9 « 0.095. 



V. THE CKM MIXING MATRIX 
HIERARCHIES (NEGLECTING CP- VIOLATION) 

A theory of flavor must explain the fermion mass hier- 
archies and the measured flavor mixing. Therefore it is 
important for the reconstruction of the Yukawa matrices 
to study if it is possible to express the absolute values of 
the CKM matrix elements as simple functions of the two 
basic flavor parameters, 9 and A. We include for com- 
pleteness a compilation of the latest extractions of the 



elements of the CKM mixing matrix, 



V, 



exp 

CKM 



0.9739 ± 0.0005 

0.2224 ±0.0035 

< 0.005 



Here to obtain \V U , 



0.2224 ±0.0036 0.00357 ± 0.00031 
0.9740 ± 0.0008 0.0415 ± 0.0008 
0.0405 ± 0.0035 0.99915 ± 0.00015 

(52) 
d \ two measurements, from superal- 
lowed Fermi transitions and nuclear beta decay, have 
been combined, as in page 36 of the 2002 CERN Work- 
shop on the CKM matrix 0. The value used for \V US \ 
was calculated in page 37 of the same reference by re- 
quiring unitarity. For \V u b\ and \V c b\ we use the latest 
extractions from B physics, as in page 6 of the same ref- 
erence. For the rest of CKM elements we use the 2002 
PDG compilation values [8( • We fit each of the measured 
absolute values of the CKM matrix elements to simple 
functions of products of integer powers of the fundamen- 
tal parameters A and 9. We use the numerical values for 
9 and A as determined from the charged lepton sector, 
i.e. 9 w 0.095 and A w 0.21, and look for correlations of 
the form r 9 p X q , where p and q are integer numbers and 
r is a low integer or rational number. We obtain as the 
best fits to the measured CKM elements the functions, 



\V US 
Kb 

*us 

Kb 

Kb 



0.2224 ±0.0036 « A 
0.187 ±0.006 «20 

0.086 ±0.009 



A 



(53) 
(54) 

(55) 



We note that the ratio |Kb/Kb| has an important uncer- 
tainly which, in principle, would allow us to fit it at la to 
both terms 9 or A/ 2. The large experimental uncertainty 
in the entry |Kd| does not allow us to implement a fit 
to 9 and A. Therefore we obtain, ignoring CP-violating 
phases, the following structure for the CKM matrix as a 



function of the parameters 9 and A, 

1 

Vckm(A,#) 



A 2 /2 -A 


a 


A 1-b 2 


-29X 


c 29X1- 


- 29 2 X 2 



(56) 



where a, b and c must be considered unknown functions 
of 9 and A which can be calculated requiring the matrix 
Vckm(A,#) to be unitary. This determines a system of 
three equations which can be solved requiring unitarity 
to order C(A 3 ), 

V ud V cd + V U8 V es +V ub V cb = 0(A 4 ), (57) 

V ud V td + V us V ts + V ub V tb = 0(A 4 ), (58) 

V td V cd + V ts V cs + V th V cb = 0(X 4 ). (59) 

We obtain a = c = 9X 2 and b = X 2 /2 + 29 2 X 2 . Therefore 
we obtain for |Vckm(A, 9)\, 



1 - A 2 /2 A 

A l-A 2 (l + 46» 2 )/2 
<9A 2 29X 1 



9X 2 
29X 
- 29 2 X 2 



(60) 



We note that if we had chosen the correlation |V^b/V^,| w 
9 instead of |14f,/T4i,| w A/2 the CKM matrix could not 
meet the unitarity requirement. 

We also note that one of the most interesting character- 
istics of the Yukawa matrices proposed in this section is 
that they account for the quark mass ratios and CKM ele- 
ments quantitatively with only two real parameters. It is 
known that a general unitary matrix can be parametrized 
using three mixing angles and a complex phase. Our re- 
sults indicate that only two real parameters seem to be 
necessary to account quite well for the absolute values 
of the CKM elements and quark mass ratios. The CKM 
reconstructed matrix in Eq. | 60 | can be expressed, using 
the standard PDG notation |8| as the product of three 
rotation matrices, each around a different axis, 



Vckm| ~ R (&i2)R 



?23/ 



)R 



13/ 



(61) 



where 9 12 = A, 9 23 = 9X and 9 13 = 9X 2 . Here R 12 (9 12 ) is 
a rotation of angle A in the first-second generation plane, 



R i2 (X) 



l-4f -A 
A l-£ 







I) 



(62) 



and R (023 ) and R (#13) are a rotations of angle 9X and 
9X 2 around the second-third and first-third generation 
planes respectively, 



R zi {29X) 



and, 



R 13 (9X 2 ) 



1 

1-29 2 X 2 -29X 
29X 1 - 26> 2 A 2 



1-^0 -9X 2 
1 
9X 2 01-^ 



(63) 



(64) 



We have seen that present experimental data allow us 
to express the fermion mass hierarchies and the absolute 
values of the CKM matrix elements as a function of two 
basic parameters 9 and A. Throughout this section we 
have ignored the presence of a CP- violating phase, which 
is required experimentally. The presence of CP-violating 
phases could affect seriously the relations between the 
absolute values of some of the CKM elements and the 
parameters A and 9 as proposed in this section. We will 
see in Sec. I VIII that CP-violating phases can be included 
in the previous analysis giving an excellent fit to the data. 
In the next section we will show that to leading order 
in A there is a unique set of Yukawa matrices that can 
reproduce these hierarchies. 



VI. QUARK YUKAWA MATRICES 
(NEGLECTING CP- VIOLATION) 

In this section we propose a particular reconstruc- 
tion of the quark Yukawa matrices that can explain the 
quark mass and mixing hierarchies found in previous sec- 
tions. We restrict our discussion to symmetric mass ma- 
trices. We will also assume that correct CKM elements 
or masses do not arise as approximate cancellations re- 
quiring the tuning of different Yukawa matrix elements. 
It is convenient to define the 3x3 normalized fermion 
mass matrices as, 



1 
m D = — m_D, 
m b 



1 
my = -mp. 
m t 



Here mo.u are the quark mass matrices and fn b and frit 
are normalized bottom and top quark masses, which are 
defined as the ratio of bottom and top quark running 
masses over the largest eigenvalue of the respective nor- 
malized matrix. We have the freedom to choose the entry 
(33) in the normalized mass matrices equal to 1, which 
correspond to the heaviest eigenvalue. Although not di- 
agonal in the gauge basis the matrix iri£) can be brought 
to diagonal form in the mass basis by a biunitary diag- 
onalization, {Vf^vaoVji = (m d ,m s ,rn b ). Analogously 
the up- type quark mass matrix, m^, can be brought to 
diagonal form in the mass basis by a biunitary diagonal- 
ization, (V^mc/V^ = (m u ,m c ,m t ). The CKM mixing 
matrix is defined by Vckm = V£ V*. 

First we note that it is not possible to generate the 
Cabibbo angle, A ~ \V US \, from the mixing between 
the first and second generations in the up-type quark 
sector. The normalized charm quark mass is approx- 
imately 9X 2 , therefore to generate the correct magni- 
tude for \V US \ from |V"| 21 , |V£| 21 « (m C /)i 2 /(mj/) 2 2, we 
would need (rhu) 12 « 9X 3 but if it this were the case 
the normalized up mass would be too heavy, m u /mt ~ 
(rhu) 12 /(™tr) 22 ~ 9X 4 , which is in disagreement with 
the measured up-top mass hierarchy, m u /m t ~ 9X e . Ob- 
taining the correct up-quark mass give us a bound on 
the entries of the upper-left submatrix of the normalized 



up-type quark matrix, which must look like, 



™ AU— c 

m U) 2 x2 



< 0(9X 6 ) < 0(6»A 4 ) 

< 0(6>A 4 ) ex 2 



(65) 



Therefore the Cabibbo angle must arise from first-second 
generation mixing in the down-type quark sector. We 
pointed out in the previous sections that the measured 
hierarchies in the down-type quark sector can be written 



as, 



m d 
m b 



9X J 



<< ^~g\ < |Ks|~A. 
m b 



(66) 



We note that assuming (nio)i2 = OX 2 and (mc)22 = 
9X « m s /m b , we obtain (m£>)i 2 /(rn.o)22 = 9X 3 w 
rrid/mb and (hid) 12/ (111.0)22 = ^- This is consistent 
with a down quark mass mainly generated from the first- 
second generation mixing as first pointed out by [3j, *-e., 



(in D 



x d—s 



ex 2 
ex 2 ex 



(67) 



since m 



case one correctly obtains |Vl| 12 



this 

(mi))i2/ '(in d) '22 =A» \V US \- Furthermore, as shown in 
the previous section, the measured hierarchies between 
the CKM elements can be written as, 



\V uh \ « OX 2 < \V cb \ « 29X ~ {m s /m b ) w 9X < \V US \ 



* A. 

(68) 
We note that if we additionally assume that (1x1.0)23 = 
29 X then we obtain the correct \V cb \, since |VfL„ = 
(m£)) 2 3/(m£i)33 = 29X. One can wonder if it would be 
possible to fully generate the CKM matrix entry \V ub \ 
from the (23)-(12) mixing in the down-type quark sector, 
or in other words, if we can assume that the normalized 
down-type quark matrix has a zero in the entry (13) and 
(31). A calculation of the diagonalization matrices show 
us that this possibility is not viable. It this were the 
case |V^| would be given by, \V[ I « 9 2 X 3 which is two 
orders of magnitude below the measured value for |V c b|, 
|V c b| ~ 9X 2 . Therefore we need to generate |V c b| directly 
from a non-zero (m£>)i3 entry, 



m D 



9X 2 9X 2 
9X 2 9X 29X 
9X 2 29X 1 



(69) 



We must note here that the possibility to fully generate 
\V ub \ in the up-sector is not viable, if that were the case 
we would generate a up-quark mass one order of magni- 
tude too heavy. This ihu matrix would predict success- 
fully all the elements of the CKM mixing matrix, assum- 
ing the mixing in the up-type sector does not affect the 
leading order predictions in A. This can be observed in 
the expressions for the diagonalization matrix, V[ , given 
by, 



(70) 





\ 2 
1 — — 

1 2 


X -9X 2 


vt = 


-A 1- 


-±L(l + 49 2 ) -29X 




-9X 2 


29X 1 - 29 2 X 2 



Using this we obtain for the mass eigenvalues to leading 
order in A, 



(V^ihnVi 



(9X 3 7 9X,1 + 9 2 X 2 ) 



(71) 



which is in perfect agreement with the reconstructed 
quark mass hierarchies in Eq. 1471 In the previous rea- 
soning we assumed that the possible flavor mixing in the 
up-type quark sector does not affect to leading order the 
predictions for the CKM matrix which are generated in 
the down-type sector. If this were the case, there are two 
simple solutions which allow us to generate the correct 
up-quark mass: directly from an entry (11), 



VClu = 



or from the first-second generation mixing, 



9X 6 








9X 2 





1 



(72) 



m v = 





9X 4 





6>A 4 





9X 2 








1 



(73) 



Since both possibilities make the same predictions for 
quark mass ratios and CKM elements both are equiva- 
lent by a rotation of the quark fields. It could be possi- 
ble to generalize the previous solution to a solution that 
generates part of | V cb \ from flavor mixing between second 
and third generations in the up-type quark sector. Let 
us assume that, 



m D 



and, 



my 



9X 2 9X 2 

9X 2 9X (2 - e)9X 
9X 2 {2-e)9X 1 



9X 6 

9X 2 -e9X 
-e9X 1 



(74) 



(75) 



In this case the diagonalization matrices are respectively, 





\2 
1 - — 

2 


-X -9X 2 


d _ 


A 1- 

(e-l)9X 2 


-^-(l + if9 2 ) r]9X 
r]9X 1 - \i] 2 9 2 X 2 



where r\ = (2 — e) and, 



VI 



l-e 2 «- 
e9X 



2 _\i 
2 





-e9X 



1 -e 



2 2 \ z 



(76) 



(77) 



These two solutions; the one given by Eos. 1691 and 1721731 
or the one given by Eg s . ITU and 1751 arc indistinguishable 



in their predictions for quark mass ratios and CKM ele- 
ments to first order in powers of A. Both reproduce the 
correct form for the CKM matrix in Eq. EH1 Therefore 
they are equivalent and can be related by a rotation of 
the quark fields. We note that from their diagonalization 
we obtain the correct empirical expressions for A and 9 as 
a function of dimcnsionless quark mass ratios (see Eas.1181 
fclT7Jl. To first order, 



m d 
m s 



m 2 b m d 



in, u 
m, 



m 2 m t 



(78) 
(79) 



Nevertheless, there are in principle other solutions that 
are not equivalent to the family of solutions here pro- 
posed and make the same predictions to leading order. 
These could be differentiated in their precision predic- 
tions for mass ratios and CKM elements when including 
higher orders in powers of A. 



general hermitian case we can introduce complex phases 
in the form, 



m D 







k 2 e^ex 2 
e-^ex 2 ex 2e^ex 

e'^OX 2 2e-^ 3 6X 1 



(80) 



By a redefinition of the phases of the quark fields this 
can be simplified to, 



m D = 



o ex 2 c^ex 2 
ex 2 ex 26X 

e^9X 2 26X 1 



(81) 



where 7 = — (V^ - V'i~'03)- If this were the case we obtain 
for the CKM matrix, Vckm = V^, to leading order in A, 



(82) 



\ 2 

1 — — 

1 2 


x -e-^ex 2 ' 


-A 1- 
(e** - 2)6X 2 


-^(l+46> 2 ) -29X 
20X 1 - 29 2 X 2 



VII. INTRODUCING CP- VIOLATION 



The angle 7 introduced this way coincides with the stan- 
dard definition, 




n/3 7t/2 



FIG. 3: Relation between CP-violating phases (3 and 7 as 
predicted by the quark Yukawa matrices in Eqs. \7<!\ and \81\ 
The horizontal hatched strip corresponds to the measurement 



sin(2/3) c 



0.78 ± 0.08. The vertical hatched strip corre- 



sponds to the la global fit for the angle 7, 7fi t = 61° ± 11°. 
The solid curve corresponds to the leading order relation be- 
tween /3 and 7 given by Eq. I #51 The dotted curved corresponds 
to the next to leading order relation given by Eq. \1U4\ 

We have seen in the previous section that a set of 
two parameter Yukawa matrices of the form given in 
Eqs. EH1 and [721 represents a family of solutions, in the 
basis where the up-type Yukawa matrix is diagonal, that 
can account for the quark mass ratios and the absolute 
values of the CKM matrix elements. It is possible to in- 
troduce complex phases in this picture to account for the 
measured CP-violation without spoiling these successful 
predictions. In order to do so we promote the real sym- 
metric matrix in Eq. 1691 to be hermitian. In the most 



7 = Arg 



V u dV* b 



VcdV c * b 
Furthermore the angle (3 is given by, 



j3 = Arg 



v td v t 



lb 



(83) 



(84) 



The angle a can be obtained from the relation, a+/3+7 = 
180°. The hermitian matrix of the form given by Ea.1811 
predicts a simple relation between the angles f3 and 7, 
which to leading order is, 



13 = Arg [2 



(85) 



The angle (3 can been determined with 7% accuracy using 
the experimental determination of sin 2/3, sin 2(3 — 0.78 ± 
0.08. We obtain f3, 



CXp 



25.7° ± 3.5°. This indicates 
that the value of (3 that nature has chosen is close to 
the maximum value that (3 can reach in the case under 
consideration, /3 max = 30°, which appears for 7 = 60°, 



d ^ n fin" 

— = — > 7 = 60 

d7 



a = 90°. 



(86) 



We note that this does not correspond with the case of 
maximal CP-violation. If we compute the determinant 

of Jarlskog matrix, C = fnr/m} r , mnfhn , we obtain, 



det C = 2J9 i X 6 



0{X X 



(87) 



where J is the Jarlskog parameter, the invariant measure 
of CP-violation, which in our case is given by, 



l 7 = 2sin(7)6> 2 A 4 (l + 0(6>A)). 



10 



We note that the maximal CP-violation case corresponds 
to 7 = ir/2 and using the Eq. [55] this corresponds to, 



dj 



= 0^7 = 90°^/3 = 26.6°. 



(89) 



Although there is a simple relation, see Eq. |SSJ between 
(3 and 7 the angle 7, as can be seen in figurc|3 cannot be 
determined with good precision from that relation and 
the experimentally determined value of [3. We obtain, 
7thco = 65° ± 38°, which is in agreement with the 2004 
winter global fit of the CKM elements, obtained using 
the results of program CKMFitter [ij , 



7cxp = 61° ± ll c 



(90) 



Alternatively we can use this experimental value of 7 to 
predict (3 from the Eq. |SSJ as can be seen on Fig. [3] We 
obtain to leading order, 



fl 



! thco 



= 29.4° ±0.05° 



(91) 



which corresponds to sin(2/?) t hco = 0.855 ± 0.001. The 
Jarlskog parameter is determined experimentally to be 
J = (3.0 ± 0.3) x 10~ 5 . The use of the Jarlskog param- 
eter does not allow us to extract 7 with better precision 
because of the uncertainties in the determination of A 
and 9. On the other hand our expression for J predicts 
an interesting relation between J and the quark masses, 



J 



2^sin( 7 )(l 



0{9\)) 



(92) 



Finally, we note that there two non-trivial characteristics 
in the relation between (3 and 7, Eq. [B3J as predicted by 
the CKM matrix in Eq. |H3 first there is no dependence 
on A or 9 to first order and second and most important 
the relation agrees with the experimental measurements 
of (3 and 7. 



VIII. PREDICTIONS FOR MASSES AND 
MIXINGS 

In this section we want to show that the simple three 
parametric Yukawa matrices proposed in Sec. lVIll can fit 
all the experimental data on the fermion spectra with 
precision. We study in this section precision predictions 
for the lighter quark masses and CKM elements. Let us 
assume that the m^ Yukawa matrix is given by Eq. OH] 
while the rhu matrix is given by Eq. [721 If we postulate 
these up- and down-type quark Yukawa matrices we can 
express the CKM elements and the quark mass ratios as 
a function of 9, A and 7. The down- type quark mass 
ratios can be calculated to the next to leading order in A 
from the diagonalization of Eg 1811 



m d 

m s 

m 



A 2 (l-(9A(4c 7 -9)), 



— = 9\{l-A9\ + Y). 



in 1, 



(93) 
(94) 



experimental input parameters 


|K.| 


0.2225 ±0.0035 


m c /m t 


(3.7 ±0.4) x 10" 3 


7 


61° ±11° 


pole 

m t 


174.3 ± 5.1 GeV 


m *( m l)ls 


4.2 ±0.1 GeV 


mv° lc 


1.7769 ± 0.0003 GeV 


predictions 


A 


0.211 ±0.007 


6 


0.083 ±0.014 


sin(2/3) 


0.824 ± 0.004 


\v ud \ 


0.975 ± 0.002 


\v ub \ 


0.0037 ± 0.0009 


\v„\ 


0.9771 ± 0.0017 


\Vcb\ 


0.035 ± 0.007 


|Vw| 


0.007 ± 0.002 


\Vu\ 


0.035 ± 0.007 


\Va\ 


0.9993 ± 0.0002 


m u (2 GeV)^ 


2.1 ±0.9 MeV 


m d (2 GeV)^ 


4.2 ± 1.4 MeV 


m s (2 GeV)^ 


84 ± 19 MeV 


m p e olc 


0.49 ±0.13 MeV 


< olC 


92 ± 17 MeV 



TABLE II: Low energy measured values of \V US \, the ratio of 
running masses m c /mt and the 1 a global fit of the phase 7 
are used to determine A and 9 from Eas. l95l fc l96l Then these 
are used to predict the three lighter quark masses, sin(2/3), 
the rest of CKM matrix elements and the muon and electron 
masses from the Yukawa matrices in Eqs. 1721 . 1811 and 11061 



Here c 7 = cos(7). The up- type quark mass ratios are 
given by, 



m c 



A 4 



9X 2 



171 1 



(95) 



while the absolute values of CKM matrix elements to the 
next to leading order in A are given by, 



\V US 

\v ud 
\v ub 
\v cs 
\v cb 
\v td 

\Vts 

\v tb 



= A-2(c 7 -2)0A 2 ±O(A 3 ), 
^A 2 ± 2(c 7 - 2)6»A 3 , 



1 



and |Vc| = |V^| 
order in A by, 



(96) 

(97) 
(98) 

2), (99) 
(100) 
(101) 
(102) 
(103) 

while (3 and 7 are related to the next 



= 6»A 2 ± 2c 7 6» 2 A 3 , 

= l-iA 2 (l±46» 2 )±26»A 3 (c. 

= 20A(1 + 0A), 

= (5-4c 7 )3 (6»A 2 ± 46> 2 A 3 ) , 

= 26»A ± 26> 2 A 2 ± (c 7 - 1)6»A 3 

= l-2\ 2 9 2 



40 3 A 3 . 



S3 = Arg [(2 - e-*T) (l ± 9\(\ - 2e^))] 



(104) 
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We use as an input three parameters determined experi- 
mentally: | V us | , m c /m t and the angle 7, whose numerical 
values can be read in table [H] We determine A, 9 solv- 
ing numerically the system of Eqs. ESI and EE1 Finally 
we use these together with the measured third genera- 
tion fermion masses to predict the masses of the lighter 
quarks, the rest of the CKM elements and sin(2/3), the 
results are shown in table[H] It is remarkable that all the 
next to leading order predictions agree with the respec- 
tive measured values. We note that the value predicted 
for sin(2/3) is slightly lower than the value predicted to 
leading order in the previous section. 

It is worth to note that a theory of flavor must predict 
succesfully the so-called Q factor. This is a combination 
of quark masses which has been determined experimen- 
tally from pseudoscalar meson masses to a 3.5 % accu- 
racy. It is defined by, 



If this were the case the charged lepton mass ratios could 
be calculated by a biunitary diagonalization, 



Q 



22.7 ±0.08. 



(105) 



In our case using the central values for 9, m tl mi> and 7 in 
tableGUand A = 0.211 ± 0.007 we obtain Q = 23.5 ± 0.80 
which agrees at la with the experimental result. For 
A = 0.218 we obtain the central value Q = 12.1 . We note 
that it is not convenient to use the measured value of Q 
to determine one of the basic parameters, instead of \y us \ 
or m c /m t , because Q contains a implicit dependence on 
the uncertainity in the top and bottom quark masses. 

To sum up, there are six input parameters in the up, 
down and charged lepton Yukawa matrices: 9, A, 7, 
mb(rrib), m\° c and m^° lc . We can choose two observables 
to determine 9 and A. In table |n] we chose |V^ S | and the 
charm/top quark mass ratio. Therefore we obtain the 
following true predictions: two quark mixing angles, a 
prediction for the CP-phase /3, the up, down and strange 
quark masses plus the electron and muon masses. A total 
of 8 true succesfull predictions as can be seen in table II. 



IX. CHARGED LEPTON SECTOR SPECTRA 
AND THE GEORGI-JARLSKOG FACTOR 

We pointed out in Sec. |H] see Eqs. El an d El that 
there are empirical relations that connect the charged 
lepton and the quark masses. In this section we argue 
that there is already a simple explanation for these rela- 
tions, the Georgi-Jarlskog factor. Let us assume that the 
normalized Yukawa matrix for the charged lepton sector 
is given by, 



9X 2 e- l ~<9\ 2 
m L = j 9X 2 3 9X 29X 
e^9X 2 29X 1 



(106) 



TO,, 



1 



A 2 (l-4(9A(cos7 



17 
12 



))• 



^ = 3*A(l-^A + in 



(107) 
(108) 



which would explain the observed empirical formulas. A 
matrix like Eq. 11061 especially the relation (01^)22 = 
3(m£>)22, could be understood in the context of grand 
unified models. For instance, by embedding the quark 
and lepton fields in the representations 5 and 10 of SU(5) 
and assuming a non- minimal Higgs structure in the uni- 
fied theory [5J such that the field that couples to the 
matter fields generating the (111^)22 entry transforms un- 
der the representation 45 of SU(5). The corresponding 
Clebsch-Gordan factors could generate the factor 3 in 
the entry (22) of the charged lepton Yukawa matrix. We 
must emphasize that even tough the original GUT model 
by Georgi and Jarslkog is ruled out the Georgi-Jarslkog 
factor, or in other words the 45 Higgs, has been used 
by many models especially supersymmetric GUT models 
which are not ruled out by current data. This may in- 
dicate that the empirical relations support a mechanism 
which can be implemented in many GUT models but it 
does not support a particular GUT model. 

Using the measured charged lepton mass ratios and the 
lcr global fit value of 7 we can determine the parameters 
A and 9 in the charged lepton sector from Eqs. 11071 and 
11081 Then using the top and bottom quark masses and 
the results of the previous section we can predict the four 
lighter quark masses, the CKM elements and sin(2/3). 
The results are presented in table Hill All the predictions 
are very close to the respective measured values, which 
is consistent with the numerical results in the previous 
section. 

Alternatively one can use the values of A and 9 as deter- 
mined from the quark sector together with the measured 
tau lepton mass to predict the electron and muon masses 
from Eqs. 11071 and 11081 The results, which are shown in 
table [HJ are consistent with the measured electron and 
muon physical masses. 



X. PERSPECTIVES AND CONCLUSIONS 

In this section we include some considerations regard- 
ing the possible characteristics of an underlying theory of 
flavor that is able to make sense of the previous results. 
To recapitulate, 

1. There are two empirical formulas that connect the 
six fermion mass ratios and the CKM elements, 

2. A simple three-parametric set of Yukawa matrices 
for the quark and charged lepton sectors can gen- 
erate these relations naturally and account for the 
measured fermion mass ratios and CKM elements, 
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input parameters 


m e /m,j. 


(4.73711 ±0.00007) x 10~ 3 


m^/rur 


(5.882 ±0.001) x 10~ 2 


7 


61° ±11° 


K olc 


174.3 ±5.1 GeV 


m b(mb)jM 


4.2 ±0.1 GeV 


predictions 


A 


0.199 ± 0.001 


e 


0.100 ±0.001 


sin(2/3) 


0.820 ± 0.005 


\v ud \ 


0.9778 ± 0.0005 


\v us \ 


0.211 ±0.003 


\v ub \ 


0.0040 ± 0.0001 


\Vcd\ 


0.211 ±0.003 


\v„\ 


0.9794 ± 0.0002 


\Vob\ 


0.0405 ± 0.0006 


\Vtd\ 


0.0075 ± 0.0009 


\Vts\ 


0.0401 ± 0.0007 


\Vtb\ 


0.99917 ± 0.00002 


m 4™c)is 


1.34 ± 0.08 GeV 


m u (2 GeV)^ 


1.85 ± 0.17 MeV 


m d (2 GeV)^ 


4.2 ± 0.3 MeV 


m s (2 GeVfe 


94 ± 5 MeV 



TABLE III: Low energy dimensionless ratios of measured 
charged lepton running masses plus the 1 a global fit value of 
the phase 7 are used to determine A and 8 from Eas. 11071 and 
11081 Then these are used to predict the four lighter quark 
masses, sin(2/3) and the CKM matrix elements from Eas. 1721 
l8Tlandrt04l 



3. The simplest explanation of the charged lepton hi- 
erarchies requires the use of grand unification to 
account for the factor 3 in entry (22) of the charged 
lepton Yukawa matrix, 

4. The proposed empirical formulas work perfectly at 
low energies but seem to get spoiled when extrapo- 
lated to very high energies, of the order of the GUT 
scale M G « 10 16 GeV. 

Additionally, the proposed Yukawa matrices have the fol- 
lowing characteristics, 

1. All the entries except entry (33) are proportional 
to a common parameter, 9, which is approximately 
9 « 0.1 

2. The generation of the correct fermion mass hierar- 
chies requires the introduction of different powers 
of A, the second flavor parameter, which is approx- 
imately A ?s 0.21 

3. The CP violating phases (3 and 7 are related by a 
simple formula, which predicts that the maximum 
value of f3 that can be reached is close to the mea- 
sured value 



We note that any theory of flavor that can generate the 
simple set of matrices proposed in this paper (or an al- 
ternative set of matrices equivalent to leading order in 
A) would automatically fit the experimental data. The 
generation of hierarchies in the Yukawa matrices, like the 
hierarchies generated by polynomial matrices in powers 
of A, is relatively easy to implement through the break- 
ing of a flavor symmetry, assuming that the vevs of the 
flavor breaking fields have a certain hierarchical struc- 
ture. There are two characteristics I want to highlight: 
the presence of a common parameter in all the entries 
of the Yukawa matrices except entry (33), and the fact 
that the empirical relations seem to get spoiled when ex- 
trapolated to very high energies. A possible theory to 
explain the presence of a common parameter in all the 
entries of the Yukawa matrices except entry (33) is the 
radiative generation of Yukawa couplings. We note that 
the parameter 9 has the right size to be a loop factor. 
If this is so, the Yukawa couplings must be generated 
at a scale not very high; otherwise our mass relations 
would be spoiled, as was pointed out above. To generate 
Yukawa couplings radiatively, one has to postulate the 
existence of additional fields belonging to two different 
sectors: the flavor breaking sector and the flavor messen- 
ger sector. The messenger sector fields would transmit 
flavor violation from the flavor breaking sector to the 
matter sector, generating Yukawa couplings radiatively. 
One more piece of the puzzle is the factor 3 in the charged 
lepton Yukawa matrix; the simplest explanation of this 
factor requires grand unification. 

Interestingly, supersymmetric models can reconcile 
the generation of the Yukawa couplings at low energy 
with grand unification 10]. It is known that grand 
unification in the context of supersymmetric models 
can successfully predict the weak mixing angle if the 
unification scale is around 10 16 GeV. This provides us 
with a consistent scenario where we can generate the 
Georgi-Jarlskog factor. On the other hand, the presence 
of soft supersymmetry breaking terms allows for the 
radiative generation of quark and charged lepton masses 
through sfermion-gaugino loops. The gaugino mass 
provides the violation of fermionic chirality required by a 
fermion mass, while the soft breaking terms provide the 
violation of chiral flavor symmetry [llj . In this case the 
superpartners of the Standard Model matter fields would 
be the flavor messengers. A supersymmetric model that 
implements the low energy radiative generation of 
Yukawa couplings has been proposed recently. This was 
achieved by postulating a U(2) horizontal symmetry 
[l2( that is broken by a set of supersymmetry breaking 
fields [T3. The model can also overcome the present 
constraints on sup ersy mmetric contributions to flavor 
changing processes 



crsy 

m 



It is known that 13 out of the 18 parameters of the 
Standard Model belong to the flavor sector: 9 fermion 
masses, 3 mixing angles and 1 CP-violating phase. We 
have shown in this paper that there are regularities un- 
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dcrlying the measured fermion masses that allow us to 
connect them through two simple empirical formulas. 
This implies a reduction in the number of fundamen- 
tal parameters in the underlying theory of flavor from 13 
to 6. We have proposed a simple set of three-parameter 
Yukawa matrices, with two real parameters and a com- 
plex phase, that can precisely account for these mass re- 
lations and give us a simpler parametrization of the CKM 
matrix The proposed set of Yukawa matrices may make 
the features of the underlying theory of flavor more ap- 
parent. Any theory of flavor that is able to generate this 
set of matrices would automatically fit the experimen- 
tal data. Furthermore, the proposed Yukawa matrices 
predict a simple and succesfull relation between the SM 
CP-violating phases. We have also pointed out that the 
empirical mass formulas between the quark and charged 
lepton masses find their simplest explanation in the con- 
text of grand unified theories. There is hope that our 
knowledge of the ligther quark masses is going to improve 
considerably in the near future by the use of lattice QCD 
methods. These empirical relations, if confirmed, could 
be a guiding light in the search for the underlying theory 
of flavor. 



Appendix 
A. Running Lepton masses 

To compute the running charged lepton masses, I use 
well known expressions, included here for completeness. 
The physical charged lepton masses are related to the 
MS running lepton masses, TOi(/x)jjg = mi(fi), through 
the relation, 



where the self-energy correction dominates the thresh- 
old. We will use s 2 w {m z )y^ = 0.23113(15), a{m z )zL = 
127.934 ± 0.027, m w = 80.423 ± 0.039 GeV and m z = 
91.1876 ±0.0021 GeV. We obtain, 

m e {m z ) = 0.487304 ± 0.000005 MeV, (116) 
"V(m z ) = 102.8695 ±0.0005 MeV, (117) 

m T (m z ) = 1748.87 ±0.30 MeV. (118) 

We note that at Q = m z the larger uncertainty in the 
running masses comes from the uncertainty in a(m z ). 
These running masses were used in Sec. [H] to search 
for correlations in higher-order dimensionless ratios of 
charged lepton masses. 



B. Running quark masses 

To calculate the dimensionless ratios of running quark 
masses we must renormalize the quark masses to a com- 
mon scale. For completeness we include in this section a 
brief explanation of the methods used to calculate run- 
ning quark masses and an update of previous numerical 
results [Tj]. Different quark masses are usually given at 
different renormalization scales. For the top quark our 
starting point is the pole mass. We use the CDF/DO 
working group average , 



m, 



174.3 ±5.1 GeV. 



(119) 



For the bottom and charm quarks we start with the run- 
ning masses, rnb(m h)jr g and TO c (m c )j^g as extracted from 
sum rules in Refs. [l£| & [16| respectively. The averaged 
values are, 



mi(n) = mf olc (1 + A/ + Az + A^) (109) 

where the one-loop self-energy correction is given by, 



A,= 



a(fi) 



I- 



mi (p) 



and the Z and W boson thresholds are given by, 



A; 



Am/ = 



«(m) 

4ttc^ 

q(m) 
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28s 2 w 



\mw 



Using the measured physical masses, 

pole _ 



m 



(110) 



(111) 
(112) 



mj— = 0.510998902 ± 0.000000021 MeV, (113) 

pole _ 

lP oie _ 1776.99 ± 0.3 MeV, 



in',;- = 105.6583568 ±0.0000052 MeV, (114) 



(115) 



we can calculate the running masses at a common scale. 
We choose to evaluate the running masses at /i = m z 



m b (m b ) m = 4.2 ±0.1 GeV, (120) 

m c (m c ) m = 1.28 ±0.09 GeV. (121) 

This value of the charm quark mass is compatible with 
recent lattice calculations, m c (m c )^L = 1.26 ±0.16 GeV 

' L ^ L/ MS 

|17j . For the lighter quarks we use the normalized MS 
values atu = 2 GeV as extracted from sum rules in 
Ref . 0, Ej . We use rescaled values @ , 

m s (2GeV)5^g = 117 ±17 MeV, (122) 

m d (2GeV) 5S g = 5.2 ±0.9 MeV, (123) 

m u (2GeV) S g = 2.9 ± 0.6 MeV. (124) 

We must add that there is a recent extraction of the 
strange quark mass by the HPQCD collaboration [201 7 us- 
ing full lattice QCD, that has extracted a central value for 
the strange quark mass lighter that the one obtained by 
sum rules and has reduced considerably the correspond- 
ing uncertainity, m s (2GeV)^L = 76±10 MeV. This value 
has not been used in the main text because it has not yet 
been confirmed by other lattice QCD collaborations. 

For simplicity and to reduce the propagation of un- 
certainties we rescale the top, bottom and charm quark 
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masses down to fj, — 2 GeV. To calculate the running top 
quark mass at fj, = 2 GeV we use the two-loop relation 
between the MS and pole quark mas ses, which is known 
through order 0(a 3 s ) [H ||l |H |H El El , 



to (m)ms 
M 

71 

H n- 

144 



= 1 + a s (n) 



>-i 



<#/*) 



3019 

' 288 






445 13 \ r 
72 36 J 

2 1 

2 + -ln2--n 

3 3 



19 n \ r , 

1 L 2 

24 12. 



S* 



(125) 



where a s (/Lt) = a s (fj,)/ir is the MS strong coupling con- 
stant, M is the on-shell mass, L = \og(M 2 / /j, 2 ) , n is 
the number of light quarks and A, A = J2i<n ( Mli> ^ s a 
small correction due to light quark mass effects [22j . 

To calculate the charm and bottom quark running 
masses at n = 2 GeV we use the analytic solution of 
the renormalization group equation in the MS scheme. 
This was originally obtained at three loops [27| and re- 
cently the four loop term has also been computed J2J| 
and found to be very small. This takes the form, 

m { \0o Po J 



2 [U> 0o 

72 To^i _ Pill + 02lo 
00 PI 



01 



a s ( M ) 2 Q26) 



where, 



70 = 1, O = (11 - §n) I, 
0i = (51 - ■ 

7i = (¥ " T»j 
/^(2857-^n-ff^)^ 

72 = (l 2 49-(^ + l|0C(3))n-^r 1 2 )i 



l 

8' 
) 16' 



(127) 



Here n is the number of light quarks, the integration 
constant to is the renormalization group invariant mass. 
We do not need to know to because, if we denote the 
right hand side as rhMsF n {^i), the running mass at scale 
/i can be calculated from a given running mass at scale 
m ( m )Ms usm S the expression, 



m (^)MS = m ( m h 



K(H) 



,MS F n (m(m))' 
In the case of four active light quarks we obtain, 

(25a s (n)Y- -"/'.. 3SD3 



(128) 



Fa{p) = 



12/25 

V 6 J 
+2.078459 a*(ji)) 



1 



3750 



a s (n) 



(129) 



To compute these we need the values of a s (mj), a s (m c ) 
and a s (n) corresponding to the experimental measure- 
ment at a s (mz)- To this end we use the three loop an- 
alytical formula for a s in the MS scheme [27]], which is 



the solution of the corresponding renormalization group 
equation, 



asM = 0Vt 



0tv 



/?ihi(t) 
01 t 

(\n(t)- 



0200 

01 



(130) 



Here t = In (/Lt 2 /A^) and A„ is the integration con- 
stant for n light quarks, to be determined from exper- 
iment. The four loop contributions to Eq. (|130f) have 
also been calculated [2£| and found to be very small. 
In practice, we use the value of a s (Mz) to first deter- 
mine A5 and a s (mj). Then we use a s (mb) to deter- 
mine A4, a s (m c ) and a s (fi — 2 GeV). Taking into 
account also the experimental uncertainty in a s (mz), 
" s (wz)ms = 0.1172 ± 0.0020, we obtain, 



A 5 


= 206 ± 26 MeV, 


(131) 


A 4 


= 277 ± 43 MeV, 


(132) 


a s (m b {m b ))^g 


= 0.218 ±0.009, 


(133) 


a, (2 GeV)^ 


= 0.286 ±0.019, 


(134) 


a s {m c {m c )) m 


= 0.354 ±0.043. 


(135) 



The uncertainty in a s (m(,(m(,)), which is four times larger 
than the uncertainty in a s (mz), comes mainly from the 
uncertainty in the determination of A5 . The uncertainty 
in the determination of A5 and A4 comes mainly from the 
uncertainties in a s (mz) and a s (mb(rrib)) respectively. 
Finally we can calculate the running quark masses at 
/i = 2 GeV. We calculate the top quark mass from 
Eq. 11251 using as an input the top pole mass in Eq. 11191 
and ais(n) as determined in Ea. 11341 We obtain, 



m t {2 GeV) 



MS 



= 298.2 ±15.4 GeV. 



(136) 



The uncertainty comes from the top pole mass uncertan- 
inty and from the the uncertainty in the determination 
of a s (fj,). Alternatively one can compute the top quark 
running mass at the top scale, m t (mt) using Eq. 11251 and 
then use formula Eq. 11281 to calculate m t {^). These two 
approaches give the same numerical results. The charm 
and bottom quark running masses are calculated from 
Eq. 11281 using as an input the running masses, m(,(mj,)j^g 
& rn c (m c )j^g, and the values of a s (mf, (to;,)), a s (m c (m c )) 
and a s (ii) determined in Eos. I133lll35l We obtain, 



to c (2 GeV)^g 
to 6 (2 GeV)^ 



1.12 ±0.13 GeV (137) 

4.91 ±0.20 GeV (138) 



Their respective uncertainties come mainly from the un- 
certainties in the theoretical extractions of Tnb(jnb)TJa 
and m c (m c )j^g in Eqs. 1120111211 These running masses 
were used together with the charged lepton running 
masses in the Sec. HH 
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